Photovoltaic effect in a gated two-dimensional electron gas in magnetic field 
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The photovoltaic effect induced by terahertz radiation in a gated two-dimensional electron gas in 
magnetic field is considered theoretically. It is assumed that the incoming radiation creates an ac 
voltage between the source and gate and that the gate length is long compared to the damping length 
of plasma waves. In the presence of pronounced Shubnikov-de Haas oscillations, an important source 
of non-linearity is the oscillating dependence of the mobility on the ac gate voltage. This results 
in a photoresponse oscillating as a function of magnetic field, which is enhanced in the vicinity of 
the cyclotron resonance, in accordance with recent experiments. Another smooth component of the 
photovoltage, unrelated to SdH oscillations, has a maximum at cyclotron resonance. 



The two-dimensional gated electron gas in a Field Ef- 
fect Transistor can be used for generation [l[ and detec- 
tion of THz radiation, and both effects were demon- 
strated experimentally [1, 0, U H 0- Concerning the 
detection, the idea is that the nonlinear properties of the 
electron fluid will lead to the rectification of the ac cur- 
rent induced in the transistor channel by the incoming 
radiation. As a result, a photoresponse in the form of 
a dc voltage between source and drain appears, which 
is proportional to the radiation intensity (photovoltaic 
effect). Obviously some asymmetry between the source 
and drain is needed to induce such a voltage. 

There may be various reasons of such asymmetry. One 
of them is the difference in the source and drain boundary 
conditions. Another one is the asymmetry in feeding the 
incoming radiation, which can be achieved either by us- 
ing a special antenna, or by an asymmetric design of the 
source and drain contacts with respect to the gate con- 
tact. Thus the radiation may predominantly create an 
ac voltage between the source and the gate. Finally, the 
asymmetry can naturally arise if a dc current is passed 
between source and drain, creating a depletion of the 
electron density on the drain side of the channel. 

The photoresponse can be either resonant, correspond- 
ing to the excitation of the descrete plasma oscillation 
modes in the channel, or non resonant, if the plasma 
oscillations are overdamped [2|. Both non-resonant [5[ 
and resonant [6|, [7| detection were demonstrated experi- 
mentally. A practically important case is that of a long 
gate, such that the plasma waves excited by the incom- 
ing radiation at the source cannot reach the drain side of 
the channel because their damping length is smaller than 
the source-drain distance. Within the hydrodynamic ap- 
proach the following result for the photoinduced voltage 
was derived for this case [2|: 



2ujt 



(i) 



where to is the radiation frequency, r is the momentum 
relaxation time, U a is the amplitude of the ac modulation 
of the gate-to-source voltage by the incoming radiation 
and Uq is the static value of the gate-to-channel voltage 
swing, U, which is related to the electron density, n, in 



the channel by the plane capacitor formula: 

en = CU. (2) 

Here, e is the elementary charge, and C is the gate-to- 
channcl capacitance per unit area. Eq. (2) is applicable 
if the scale of the variation of the potential in the channel 
is large compared to the gate-to-channel separation. 

Recently, the first experiments on the photovoltaic ef- 
fect at terahertz frequencies in a gated high mobility two- 
dimensional electron gas in a magnetic field were per- 
formed [1, Q. The main new results are: (i) the pho- 
toinduced dc drain-to-source voltage exhibits strong os- 
cillations as a function of magnetic field, similar to the 
Shubnikov-de Haas (SdH) resistance oscillations, and (ii) 
the oscillation amplitude strongly increases in the vicin- 
ity of the cyclotron resonance. 

In this Letter we consider theoretically the photo- 
voltaic effect in a gated electron gas in a magnetic field 
assuming, as in Ref. Q, that the incoming radiation 
creates an ac voltage between the source and the drain. 
Further, in accordance with the experimental conditions 
we assume that 1) the source-drain length, L, (x direc- 
tion) is greater than the plasma wave damping length, 
so that the plasma waves excited near the source do not 
reach the drain, and 2) the sample width, W, in the y 
direction, is much greater than L, see Fig. 1. The first 
assumption means that the boundary conditions at the 
drain are irrelevant and, as far as plasma waves are con- 
cerned, the sample can be considered to be infinite in 
the x direction. The second one implies a quasi-Corbino 
geometry (all variables depend on the x coordinate only). 

We explain the observed strongly oscillating photore- 
sponse as being due to the non-linearity originating from 
the oscillating dependence of the mobility on the Fermi 
energy, and hence on the ac part of the gate voltage. 

The photovoltaic effect is due to a radiation-induced 
force G driving the electron current. Without mag- 
netic field, G is obviously directed in the x direction and 
is compensated by the appearance of an electric field. 
In the presence of magnetic field the problem becomes 
more subtle, not only because in this case G has a y- 
component, but also because this radiation- induced force 
becomes non-potential: curlG ^ 0. The non-potential 
part will drive an electric current along closed loops. 
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FIG. 1: Assumed design and geometry. The THz radiation 
produces an ac voltage U a between the source and the gate 
inducing a dc source-drain voltage AC/. The gate width, W, is 
much larger than the gate length L (quasi-Corbino geometry) 



The significance of the cyclotron resonance for the pho- 
tovoltaic effect is related to the well-known dispersion 
relation for plasma waves in a magnetic field 10]. For 
gated two-dimensional electrons it reads: 



uj = + s 2 k 2 



(3) 



where to c is the cyclotron frequency, s is the plasma wave 
velocity, and k is the wavevector. Thus, the plasma waves 
can propagate only if lu c < lu. In the opposite case the 
wavevector becomes imaginary, so that the plasma oscil- 
lations rapidly decay away from the source. The change 
of regime when the magnetic field is driven through its 
resonant value will manifest itself in the photoresponse. 

Following Refs. p], and other theoretical work, we 
will use the hydrodynamical approach because, like the 
Drude equation, it provides a relatively simple descrip- 
tion, compared to the full kinetic theory. However it 
should be understood that at low temperatures, at which 
the experiments [1, Q were done, this approach strictly 
speaking is not justified, because the collisions between 
electrons are strongly suppressed by the Pauli princi- 
ple. Nevertheless, the qualitative physical results de- 
rived from the kinetic equation and from the hydrody- 
namic equations are usually similar, e.g. the properties 
of plasma waves are identical in both approaches, pro- 
vided that the plasma wave velocity s is greater than the 
Fermi velocity, so that the Landau damping can be ne- 
glected [ll|. For this reason, we leave the much more 
complicated approach based on the kinetic equation for 
future studies. 

The electrons in a gated 2D channel can be described 
by the following equations: 

3~V 6 C 

— + (v-V)v = VC/+ — Bxv-7v, (4) 

at m mc 

^+div(C/v)=0, (5) 

where v is the electron drift velocity, B is the magnetic 
field along the z direction, m is the electron effective 



mass, and 7 = 1/r. The parameter 7 is an oscillating 
function of the electron concentration (or gate voltage) 
and magnetic field, which results in the SdH oscillations. 

Eq. (4) is the Euler equation, taking account of the 
Lorentz force and damping due to collisions, ft differs 
from the conventional Drude equation only by the con- 
vective term (v • V)v. Equation (5) is the continuity 
equation rewritten with the use of Eq. (2). 

The boundary condition at the source (x = 0) is: 



(7(0, t) = C7 + U a cos cot, 



(6) 



where u> is the frequency of the incoming radiation, and 
C7 a is the amplitude of the radiation-induced modulation 
of the gate-to-source voltage. For a long sample, the 
boundary condition at the drain is 



0, C7 -> C7 for x 



(J) 



We will search for the solution of Eqs. (4) and (5) as 
an expansion in powers of C7 Q : 



vi + v 2 , C7 = C7 + Ux + U 2 . 



(8) 



Here Vi and U\ are the ac components proportional to 
C7 a , which can be found by linearizing Eqs. (4, 5), v 2 
and t/ 2 are the dc components, proportional to C7^ (we 
are not interested in the second harmonic terms ~ U 2 ). 
It is convenient to introduce u = eU/m, u a — eU a /m, 
and the plasma wave velocity in the absence of magnetic 
field s = u Q /2 = (eUo/m) 1 / 2 0). 
To the first order in C7 a , we obtain: 

dv lx du\ 



dt 



+ -7^7 + VcVly + JV lx = 0, 



~dt 



U C V lx + JVly = 0, 



, dvi x 
dx 



0. 



(9) 



(10) 



(11) 



where oj c = eB/mc is the cyclotron frequency. The 
boundary conditions follow from Eqs. (6, 7): ui(0, t) = 
u a cos(wt) and ui(oo, t) = 0, Vi(oo, t) = 0. 

Searching for the solutions ~ exp(ikx — iut), we obtain 
the dispersion equation for the plasma waves: 



1 + ia 



7 2 



1 



(12) 



where a = (ujt)^ 1 and f3 — uj c /uj is the magnetic field in 
units of its resonant value for a given lu. To ensure the 
boundary condition at x — > 00 the root with a positive 
imaginary part of k should be chosen. If damping is 
neglected (a = 0), this equation reduces to Eq. (3). 
The explicite expressions for U\, v± x , and v\ y are easily 
obtained from Eqs. (9-11). 
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In the second order in U a , we find 
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where the angular brackets denote the time averaging 
over the period 2it/lu. Here we have expanded the func- 
tion 7(w) to the first order in u%. The quantities 7 
and 7' = dj/du should be taken at u = uq. The 
boundary conditions for Eqs. (13-15) are: 1*2(0) = 0, 

V2x{0O) = V2y(oo) = 0. 

From Eq. (15) we derive the obvious fact that j x = 
(j x differs from the x component of the true current 
density only by a factor en). Using this, and introducing 
the y component of the current, j y , by a relation similar 
to Eq. (15), we can rewrite Eqs. (13, 14) as follows: 



Ucjy = G x (x) 



du 2 
dx 



jj y = G y {x), (16) 



where the additional driving force G induced by the in- 
coming radiation is given by: 



G x 



G y — 
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Uq 



-i < 
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dv\ x 
dx 



>, (17) 



7' I (u\vi v ) -(uivix) - ( 



Uq 



i'Li 



dx 



(18) 



Both G x and G y depend on x as exp(— 2k" x), where k" 
is the imaginary part of the wavevector defined by Eq. 
(12), reflecting the decay of the plasma wave intensity 
away from the source. Thus curlG ^ 0. 

One could solve Eqs. (16) to obtain the photoinduced 
voltage Au = J Q [G x — (uj c /^)Gy} dx and this would be 
the correct result for the true Corbino geometry, where 
the current j y can freely circulate around the ring. How- 
ever, we believe that this is not correct for a finite strip, 
even if W >> L, because in this case the current j y 
induced by the non-potential part of the driving force, 
G y (x), obviously must return back somewhere, forming 
closed loops [12]. How exactly this will happen, is not 
quite clear. In our model, the current loops are likely 
to close through the source contact, however in reality 
the oppositely directed y-current will probably flow in 
the ungated part of the channel adjacent to this contact. 
Anyway, since the current j y integrated over x must be 
zero (except near the extremities), we believe that the 
correct way is to integrate the first of equations (16) tak- 
ing this into account, and to ignore the second one, which 
is not applicable beyond the gated part of the channel. 
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FIG. 2: The functions f(/3) (left) and g((3) (right) describing 
respectively the smooth part and the envelope for the oscil- 
lating part of the photovoltage. The values of the parameter 
a = (lot)- 1 : 1 - 0.2, 2- 0.4, 3- 0.8 



The integration interval should be expanded to include 
the region where the current lines return backwards. 

So far, we have no rigorous proof that this idea is cor- 
rect, however we have checked that both methods give 
similar qualitative results (but differ in the exact form of 
the magnetic field dependence of the photovoltage) . 

As described above, we obtain Au = ^2(00): 



Au 



G x (x) dx. 



(19) 



Using Eqs. (17, 19) we finally calculate the dc photovolt- 
age AU — mAu/e, between drain and source induced by 
the incoming radiation: 



AU 



4 0b 



f(p) 



dry n 
dn 7 



(20) 



Here we have separated the photoreponse in a smooth 
part and in an oscillating part. The second one, propor- 
tional to dry/dn, is an oscillating function of gate voltage 
or magnetic field ~ dp xx /dn, where p xx is the longitudi- 
nal resistivity of the gated electron gas. 

Note, that even if the amplitude of the SdH oscillations 
is small, the parameter \dp xx / dn\(n / p xx ) can be large, so 
that the oscillating contribution may dominate. 

The frequency and magnetic field dependences of the 
photovoltage are described by the functions f(§) and 
g((3), which are given by the following formulas [l3j: 



f(P) 



1 



1 + F 

Va 2 + F 2 ' 



m = 



l + F 



y/a? 



F 2 



(21) 



(22) 



where F depends only on the ratio f3 — lo c /uj and the 
dimensionless parameter a — (wt) _1 : 



F 



P 2 



(23) 
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Reduced magnetic field, (3 

FIG. 3: Magnetic field dependence of the oscillating part of 
the photovoltage for a = 0.1, LOT q = 0.5 for two values of N = 
Ep/huj. The vertical scale for the lower trace is expanded 4 
times with respect to the upper trace. (5 — uj c /lo 



In the absence of magnetic field, (3 = 0, F = 1, and Eq. 
(21) reduces to Eq. (1). 

Figure 2 shows the behavior of the functions f((3) and 
g(/3) for several values of the parameter a. One can see 
that for small values of a (or large cut) the smooth part 
displays the cyclotron resonance with the unusual line- 
shape f((3) ~ [(1 - /3) 2 + a 2 ]- 1 / 2 . The envelope for the 
oscillating part, g(/3) exhibits a fast decay beyond the 
cyclotron resonance (f3 > 1), confining the oscillations of 



the photovoltage to the region f3 ~ 1. 

To display the oscillating contribution, we take the pa- 
rameter 7 in the conventional form [l4T |. which is valid 
when the SdH oscillations are small: 



7 = 7o 



X 



sinhx 



exp 



(2ttE f \ 



(24) 



where \ — 2ir 2 kT /Huj c , r q is the "quantum" relaxation 
time, and Ep is the Fermi energy, which is proportional 
to the electron concentration n, and hence to the gate 
voltage swing U. 

We introduce the parameter N = Ep/fujJ, which is the 
number of Landau levels below the Fermi level at cy- 
clotron resonance. Figure 3 presents the oscillating part 
of the photovoltage [the function — (drf/dn) (n/j)g 
for a = 0.1, x = 0.7 (corresponding to T = 4K, 
lo = 2n ■ 2.5 THz), and ujr q = 0.5, for two values of N. 

In spite of the crudeness of our model, which does not 
account for various features of the experimental situation 
(the unavoidable presence of ungated parts of the chan- 
nel, etc), our results show a good qualitative agreement 
with the recent experimental findings Q. 

In summary, we have calculated the photovoltage in- 
duced in a gated electron gas by THz radiation in the 
presence of the magnetic field. As a function of magnetic 
field, the photoresponse contains a smoothly varying part 
and an oscillating part proportional to the derivative of 
the SdH oscillations with respect to the gate voltage. The 
smooth part shows an enhancement in the vicinity of the 
cyclotron resonance. 

We appreciate numerous helpful discussions with Woj- 
ciech Knap, Nina Dyakonova, Maciej Sakovicz, Stephanc 
Boubanga-Tombet, and Sergei Rumyantsev. 
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